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ABSTRACT 


Prior  studies  of  IMPATT  diodes  have  restricted  their  consideration 
to  one-dimensional  representations  of  these  devices.  In  the  present 
report,  some  qualitatively  new  modes  of  behavior  are  brought  to  light 
from  a  study  of  a  more  realistic  multidimensional  structure.  The  pri¬ 
mary  result  of  this  study  is  the  uncovering  of  a  varied  spectrum  of  new 
avalanche  modes,  whose  character  is  discussed  in  some  detail. 

This  spectrum  is  found  to  be  comprised  of  two  different  types  of  radial 
TM  modes  which  are  referred  to  as  the  symmetric  and  the  unsymmetric 
modes.  The  totality  of  prior  results  of  one -dimensional  studies  turn 
out  to  constitute  the  lowest  order  mode  of  the  symmetric  mode  group. 
This  is  a  quasi-TEM  radial  wave  avalanche  mode,  which  is  a  bulk  mode. 
The  remaining  symmetric  modes  that  have  been  investigated  turn  out  to 
constitute  a  new  type  of  avalanche  phenomena  involving  avalanche  sur¬ 
face  wave  modes.  The  unsymmetric  modes  correspond  most  simply  to 
bulk  wave  modes,  but  differ  considerably  from  the  quasi-TEM  bulk 
mode  mentioned  above. 
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THE  MODE  SPECTRUM  OF  AVALANCHE  DIODES 


L  INTRODUCTION 

Recent  experimental  and  theoretical  investigations  have  revealed  the  existence  of  two- 
dimensional  effects  in  certain  Gunn  effect  and  other  semiconductor  structures,^'^  These  ef¬ 
fects  usually  involve  radial  variations  of  the  RF  fields  in  the  semiconductor,  in  addition  to  the 
expected  axial  variations  parallel  to  the  applied  DC  electric  field.  In  this  report,  new  micro- 
wave  frequency  modes,  with  radial  variations,  are  reported  for  IMPATT  diodes  biased  by  an 
axial  DC  electric  field. 
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Fig.  1.  Multiple  uniform  layer  model  of  avalanche 
diode.  In  the  present  work,  the  axial  lengths  of 
the  drift  zones  are  assumed  to  be  negligible. 

The  device  characteristics  for  a  uniform,  right-circular,  cylindrical,  semiconductor  pill 
of  axial  length  L  and  radius  R  (see  Fig.  1)  can  be  inferred  from  the  avalanche  phenomena  de¬ 
scribed  by  Eqs.(1)  and  (2)  for  electron  and  hole  current  continuity,  and  by  the  Maxwell  field  equa¬ 


tions,  Eqs.  (3)  and  (4). 

an/at  =  -(1/2)  V  •  J n  +  g  >  (1) 

ap/at  =  (1/2)  V  .  +  g  ,  (2) 

vxE=PoaH/at  ,  (3) 

vxH  =  j  +J  +€  aE/at  .  (4) 

n  p  '  '  ' 


Here  n  and  p  are  the  electron  and  hole  densities,  =  +qpv  and  =  — qnv  are  the  hole  and 
electron  current  densities  which  are  assumed  to  flow  at  the  same  saturated  drift  speed  |v|  par¬ 
allel  to  the  axial  direction,  and  e  are  the  permeability  and  permittivity  of  the  material,  H  is 
the  magnetic  field,  E  is  the  electric  field,  V  is  the  del  operator,  g  =  a|v|(n  +  p)  is  electron- 
hole  generation  rate,  and  a  is  the  ionization  rate  which,  for  simplicity,  has  been  assumed  to 
be  equal  for  holes  and  electrons  as  in  Refs.  5  and  6. 

To  simplify  matters,  analysis  is  restricted  to  small -signal,  time-harmonic  fields  in  the 
single  layer  case  of  distributed  avalanche,  as  in  Ref.  5,  with  =  0  at  z  -  L/2,  and  .1^  =  0  at 
z  =  — L/2.  The  influence  of  the  highly  conductive  regions  sandwiching  the  avalanching  space- 
charge  layer  is  approximated  by  using  the  boundary  conditions  appropriate  to  infinitely  conduc¬ 
tive  material,  i.e.,  the  radial  component  of  E  is  zero  at  z  =  ±L/2. 
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II.  l‘MKLDS  AND  CUHHDNTS 


It  can  be  shown  that  the  possible  field  modes  are  similar  in  some  respects  to  the  TK  and 

7 

TM  radial  waveguide  modes.  When  the  circumferential  variation  of  the  fields  and  currents  is 
assumed  to  be  negligible  (which  is  consistent  with  the  usual  device  symmetry),  these  mode  types 
are  uncoupled  from  each  other.  In  addition,  the  TE  radial  modes  are  uncoupled  to  the  avalanche 
currents  and,  hence,  are  the  same  as  for  a  "cold"  dielectric  cylinder  bounded  by  metal  caps. 
Thus,  primary  interest  is  in  the  TM  radial  modes  for  which  the  basic  form  for  the  axial  compo- 
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nent  of  E  is 
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Here  d^(x)  is  the  zero^*^  order  Bessel  function  of  the  first  kind,  T  is  the  radial  wavenumber, 
and  are  the  axial  wavenumbers,  and  A  and  B  are  constants.  Solution  of  the  boundary 
value  problem,  posed  above,  yields  a  discrete  spectrum  for  T,  K^,  and  which  in  turn  implies 
the  same  for  B.  Each  set  of  associated  values  for  T,  K^,  and  corresponds  to  a  particular 
field  and  current  distribution  (which  varies  with  frequency)  within  the  diode  that  is  referred  to, 
for  simplicity,  as  an  avalanche  mode.  Mode  is  understood  here  to  have  the  same  meaning  as, 
for  example,  used  for  empty  radial  waveguides.  The  relative  amplitudes  of  the  avalanche  modes 
would  be  determined  by  imposing  boundary  conditions  at  the  "rim"  of  the  semiconducting  pill, 
i.e.,  at  r  =  R,  in  the  same  manner  as  is  done  for  radial  waveguides. 

The  TM  fields  and  currents,  in  the  avalanche  layer  under  consideration,  have  been  shown 
to  be^ 


E^  =  J^(Tr)  2A  {cos(K^z)  +  B  cos  (K^z)} 
H  =  2iaJ6TAJ,(Tr 

(p  1' 


cos(K^z)  Bcos(K^z) 
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E  =  -2JTAJ,(Tr) 
r  1 


sin(K^z)  BK^  sin(K^z) 
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J  =  J  (Tr)  A 
n  o 


(J ^  cos(K^z) - — sin(K^z) 


+  B 


jVK  (7 

(j^  cos(K-z) - sin(K^z) 
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J  =  J  (Tr)  A 
p  o' 
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(8) 
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From  the  boundary  condition  =  0  at  z  =  L/2  (which  is  equivalent  to  =  0  at  z  =  — L/2  because 
of  the  symmetry  in  the  present  case),  Eq.  (10)  yields 


CO  cos0^  +  jVK^  sine^j 
CO  cos0^  +  jVK^  sin0^J  lc7^J 


(11) 
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where 


-2ja;V»'  I 
o  o 


‘  +  2ja)0  V  -  (K.V)^ 

•'  o  1 


(i  =  1.  3), 


©1  H  K^L/2, 


-  K^L/Z, 


O'  -  ionization  coefficient  evaluated  at  K  equal 
^  to  the  DC  bias  value,  E 

o" 

o'  =  do/dE  evaluated  at  E  =  E  , 
o'  o 

cj  -  radian  frequency, 

I  =  DC  current  density  due  to  E  . 
o  o 


(12) 


HI.  DISPERSION  RELATIONS 


The  dispersion  relation  for  the  model  waves  in  the  avalanche  diode  can  be  obtained  for 
Eqs.(8)  and  (11)  and  the  boundary  condition  E^  =  0,  at  z  -  ±\Jz,  The  resulting  dispersion  rela¬ 
tion  is  obtained  in  the  simultaneous  system  of  equations: 

F(e^)=F(e3)  ,  (13) 

where 


F(e)  E 


0  sine 

2  2 


1  -[e/e^]^  +  j/e^ 

^s  0  +  [je  sine]/0^ 


and 


0^  =  WL/2V, 


(e^  3)^  =  I  *  ^/b^  -  4C  ,  (14) 

where 

b  £  +0^-0,^,  0rp  E  TL/2, 

c  E  0^0^  -0^  [0^  +  0^ I. 
o  m  1  m  a 

2  2  2 
0  E  0  f  -  0  +  10  , , 

m  d  a  d' 

0^^  E  [L/2V]^  2o'I  V/€. 
a  ^  ^  00' 


There  turn  out  to  be  two  kinds  of  solutions  to  the  above  set  of  equations.  The  first  kind, 
which  shall  be  called  symmetric  for  reasons  to  be  discussed  shortly,  is  given  approximately  by 

T  ^j2n7r/L  ,  (15) 

n  j  / 

2  2  2 

and  Eqs.  (12)  with  Eq.  (14)  used  with  0,^,  replaced  by  — n  tt  .  These  analytical  results  have  been 
checked  with  computer  calculations  for  the  modal  index  number  n  =  0,  ±1,  ±2,  ±3  and  were  found 
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to  be  accurate  to  within  5  percent  for  the  usual  value  of  device  parameters  within  the  microwave 
range.  It  is  not  clear  at  present  how  accui*ate  these  results  are  for  n"  ±4,  However,  it  can  be 
shown  that  they  definitely  fail  as  n  goes  to  infinity  as  follows.  From  Eq,  (13),  as  0,^,  — 

0^^  — 2b — ® 

03^-c/2b-je^  . 

Substituting  these  results  into  Eq.  (13),  a  contradiction  is  obtained  because 
i/jed  ^  ^  ■*  [  1  cot 
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An  important  simplification  occurs  in  the  n.^  3  cases  if  l^j^l  »  which  is  often  valid 
in  practice.  Then, 


Mn 


4n^7r^(e^  +  je^) 

^2  .  2  2v  T  2 

(Bm  +  n  7r  )  L 
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^  2  2,  ,  2 

(©m  +  n  TT  )  L 


(17) 


IV.  SURFACE  WAVE  MODES 


The  following  may  be  observed  for  the  approximate  dispersion  relations  given  by  Eqs.  (15), 
(16),  and  (17): 


a.  Wlien  n  =  0,  then 


(18) 

Mo  ° 

(19) 

(20) 

This  mode  very  closely  approximates  a  TEAI  radial  wave  avalanche  mode  and  will  be  re¬ 
ferred  to  the  quasi-TEM  mode.  T^  niay  be  estimated  more  accurately  in  the  present  case  by 
the  result  for  an  exact  TEM  radial  wave  avalanche  mode  (which  cannot  exactly  satisfy  the  bound¬ 
ary  conditions  but  does  satisfy  them  to  a  good  approximation), 

Tq  =  [(Jq  +  jwe]  ,  (21) 

where  =  o'j(K.  =  0).  The  results  of  one-dimensional  theory^  turn  out  to  constitute  this  quasi- 
TEM  mode.  It  should  be  noted  that  this  mode  is  a  bulk  avalanche  mode  as  distinguished  from 
surface  wave  avalanche  modes  which  are  discussed  next. 


b.  The  first  higher  order  mode  is  n  =  ±  1  for  which 
T  2  3  -(Ztt/L)^  , 


(22) 
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The  qualitative  aspects  of  the  behavior  for  the  higher  order  modes,  with  n  =  ±  2,  -t  3,  are 

the  same  as  the  n  -  i  1  mode.  Here  the  radial  wavenumber  (T.)  is  of  the  same  order  of  mag- 

^  4-1 

nitude  as  K^.  which  for  typical  device  structures  is  of  the  order  of  1  0  cm  In  addition,  K. 

31  1 

is  no  longer  close  to  zero  (since  it  involves  a  factor  of  L  ),  although  it  remains  smaller  than 
and  satisfies  the  simple  relation 


+  (2K^^]2  -  +  [27r/L]2  .  (25) 

c.  The  total  current,  defined  as  Jrp  ^  ^  ^p  ^  varies  rapidly  with  axial  position,  z, 

within  the  sample  length,  L.  This  is  in  contrast  to  the  n  =  0  case  for  which  is  essentially 
independent  of  z. 

d.  Because  of  the  rapid  decrease  of  field  strength  and  current  with  radial  distance,  these 
quantities  at  a  distance  of  L  are  substantially  less  than  those  at  the  surface  of  an  avalanching 
semiconductor  pill.  Thus,  it  is  appropriate  to  characterize  the  n  >  1  avalanche  modes  as  sur¬ 
face  modes. 


V.  BULK  MODES 

In  this  section,  the  second  category  of  solutions  to  the  dispersion  relations,  Eqs.  (12)  and 
(14),  are  discussed.  These  will  be  referred  to  as  unsymmetric  modes. 

The  unsymmetric  modes  are  obtained  from  the  dispersion  relation  solution 

0^  =  03  (26) 

and  can  be  determined  in  closed  analytical  form.  The  status  of  these  solutions  is  not  completely 
clear  because  0^  -  03  implies  =  K3.  This  constitutes  a  degenerate  case  for  which  the  deriva¬ 
tion  of  the  dispersion  relations  is  not  rigorously  valid.  Those  modes  are  of  physical  interest 
because  they  represent  plausible  solutions  for  the  unsymmetric  case  where  the  saturated  elec- 


the  electron,  a  ,  does  not  equal  the  ionization  coefficient  for  the  holes,  a  )  in  the  limit  as 
n  ^  p 

V  —  V  and  a  —a  approach  zero, 
n  p  n  p 

The  solutions  are  plausible  because  of  the  continuity  normally  associated  with  physical  prob¬ 
lems,  and  thus  anticipated  in  the  present  one.  However,  a  formal  analysis  of  the  unsymmetric 
case  (i.e.,  where  ^  V^,  and/or  ^  cv^)  is  required  for  rigorous  justification.  In  present 
day  devices  ^  ^p^  small  but  not  zero  (typically  0.002  -^  0.004),  and  —01^]/ 

+  ttp]  is  small  for  Ge  and  GaAs  devices.  Thus,  these  unsymmetric  modes  may  be  of  consid¬ 
erable  physical  interest. 

2 

The  exact  unsymmetric  modal  solutions  are  determined  from  Eq.(14)  by  setting  b  —  4C 
equal  to  zero.  These  solutions  are 
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(27) 


T 


1,2, 3, 4 


±jK  1  ±  li  +  [K^  -K^]/K 
‘’a  V  ^  m  o 


a 
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Numerical  evaluations  of  the  above  modal  solutions  possess  radial  variations  that  are  relatively 
moderate,  compared  to  the  previously  mentioned  surface-wave  avalanche  modes,  and  thus  may 
be  considered  a  form  of  bulk-mode  waves  (see  Figs.  2  and  3).  However,  these  radial  variations 
are  not  negligible,  as  they  are  in  the  more  familiar  quasi-TEM  mode.  These  unsymmetric 
modes  correspond  more  closely  to  the  familar  bulk  modes  for  frequencies  in  the  lower  portion 
of  the  microwave  spectrum. 


Fig.  3.  Imaginar  art  of  the  radial 
wavenumber  Im  vs  radian  fre¬ 
quency  CJ. 


Fig.  2.  Real  part  of  the  radial 
wavenumber  Re  (T}  vs  radian 
frequency  a;. 


The  magnetic  fields  can  be  determined  from  the  description  of  E,  for  all  the  modes  dis¬ 
cussed,  from  Eq.  (3).  From  this  information  the  impedance  of  an  IMPATT  diode,  operating  in 
any  mode  or  combination  of  modes,  can  be  determined  in  a  manner  which  allows  for  the  finite 
radius  of  the  diode.  This  is  discussed  further  in  the  Appendix, 

VI.  CONCLUSION 

A  variety  of  new  modes  of  behavior  have  been  uncovered  from  a  detailed  examination  of  a 
multidimensional  avalanche  diode  structure.  The  totality  of  results  from  prior  one -dimensional 
studies  constitute  only  the  lowest  order  mode  of  one  mode  group.  The  remaining  modes  exhibit 
qualitatively  different  characteristics  from  the  one-dimensional  like  mode. 

The  potential  existence  of  these  higher  order  new  modes  might  be  anticipated  from  the  obser¬ 
vation  that  the  wavelength  of  the  fields  in  the  prior  one- dimensional  studies  is  of  the  same  order 
of  size  as  the  axial  dimension  of  the  avalanche  zone  (see  Fig.  4).  However,  it  may  be  initially 
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Fig.  4.  Axial  wavelength  of  fields  in  a  1 -micron 
avalanche  zone  in  an  avalanche  diode  determined 
from  \  =  Re  {2ir/K  }  . 

surprising  that  the  radial  variations  for  these  new  modes  are  so  rapid.  Because  the  electron 

and  hole  currents  are  assumed  to  be  solely  in  the  axial  direction  (and  interact  strongly  with  the 

axial  component  of  electric  field),  it  is  not  unexpected  that  there  are  longitudinal  variations  of 

fields  and  currents  very  different  from  those  in  the  comparable  cold  dielectric.  On  the  other 

hand,  it  might  be  guessed  that  the  radial  variations  of  fields  and  currents  should  approximate 

those  found  in  the  cold  dielectric  case.  The  fact  that  the  radial  variations  are  in  fact  rapid  over 

distances  of  only  several  microns  is  not  peculiar  to  the  avalanche  effect  and  arises  because  the 

modes  are  TM  and  not  TEM.  Waves  with  comparable  rapid  radial  variations  have  been  very 

1 

recently  discussed  by  Kino  and  Robson  for  Gunn  effect  devices,  and  by  Burke  and  Kino  for  InSb 

3 

in  the  presence  of  transverse  magnetic  fields. 

The  origin  of  the  rapid  radial  variations  for  the  radial  TM  avalanche  modes  can  be  math¬ 
ematically  traced  to  the  dispersion  relation  written  simply  as 

=  [K^  -  K^]  [1  -  f(K)]  .  (28) 

For  the  quasi -TEM  mode,  f(K)  is  practically  one  and  hence  T  is  usually  negligible.  For  the 
TM  modes,  f(K)  differs  considerably  from  unity  and  T  varies  essentially  as  K  n/1  —  f(K). 

For  conventional  IMPATT  diodes  (as  opposed  to  sheet  or  doughnut  diodes),  one  may  antic¬ 
ipate  that  in  conventional  applications  the  quasi-TEM  bulk  mode  will  dominate  all  the  surface 
wave  modes  in  device  performance.  Recent  small-signal  measurements  agree  with  this  assump¬ 
tion  and  indicate  that  either  the  unsymmetric  modes  are  not  valid  (for  reasons  discussed  earlier) 
or  do  not  couple  well  to  the  microwave  cavity  modes. 
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APPENDIX 
DEVICE  ACTIVITY 


In  generaE  for  any  device  the  electromagnetic  power  delivered  to  its  environment  is  given 

by 


for  time-harmonic,  single -frequency  fields,  where  the  integration  is  oyer  the  volume  of  the 
device.  For  the  IMPATT  device  model  considered  in  this  study,  J  has  only  a  z -component  so 
that 


-P  =  Re 


rd(^rdz 


(30) 


where  rd</>drdz  is  the  element  of  volume  in  cylindrical  coordinates. 
For  the  unsymmetric  bulk  IMPATT  modes 


-P  ^  A  Re 


r 


+  L/2  E  J* 
'  z  z 


L/2 


dz 


(31) 


where  A  =  ttR  and  the  equality  is  approximate  because  the  radial  variations  are  neglected.  It 

is  useful  to  introduce  here  the  z -component  of  the  total  RF  current  density,  J_  =  J  +  j  +  iojcE 

1  n  p  •'  z 


Then, 


^  A  Re 


+  L/2 

-L/2  ^ 


dz 


(32) 


It  is  further  useful  to  observe  that  taking  the  divergence  of  Eq.  (4)  yields 


V  .  z:  dJ^/dz  +  dJ^/dr  -  0 


(33) 


where  is  the  radial  component  of  the  total  RF  current  density  vector  Since  there  are  no 

radial  conduction  currents  in  the  assumed  model,  J  =  iwcE  ,  so  that 

r  r 


(34) 


It  follows  from  Eqs.  (3),  (4),  and  (5)  that  the  radial  dependence  of  E^  is  J^(Tr)  where  J|(x)  is  the 
first-order  Bessel  function  of  the  first  kind.  For  moderate  radial  variations,  |T|r  is  small 
compared  to  unity  when  r  is  not  "too  large."  For  that  case,  J^(Tr)  =  Tr/2  and  hence  8E^/0r  is 
independent  of  r.  Thus,  Eq.  (34)  may  be  written  as 


a(z) 


(35) 


where  a(z)  is  independent  of  r.  A  partial  integration  of  Eq.  (35)  yields 
=  b(z)  +  g(r) 


(36) 


9 


where  b(z)  =  /  a(z)  dz  and  g(r)  is  independent  of  z.  On  the  other  hand,  it  follows  from  the  def¬ 
inition  of  Jrp  and  Eqs.  (6)  through  (10)  that  for  any  mode  is  a  product  of  .J^(Tr)  and  some  func¬ 
tion  of  z.  The  only  way  these  two  results  for  the  form  of  J,p  can  be  consistent  is  if 

Jrp  =  constant  ,  (37) 

where  the  constant  is  independent  of  r  and  z.  Hence,  Eq.  (32)  becomes 

-P  S  A  |J^|2  Re{Z}  ,  (38) 

where  the  impedance  is 


Z  =  V/J^ 


(39) 


and  where  the  voltage  V  is  given  by 
p+L/2 

-V  =  \  E  dz  .  (40) 

-'-L/2  ^ 

The  preceding  argument  demonstrates  that,  if  the  radial  variations  are  sufficiently  small, 
the  power  output  from  radial  TM  modes  is  given  by  the  conventional  result  for  which  the  conven¬ 
tional  derivation  implicitly  considers  only  a  radial  TEM  mode.  Significant  errors  will  arise, 
however,  where  |T|  r  begins  to  be  comaprable  to  unity.  Since  E^  and  vary  as  J^(Tr),  a 
correction  to  Eq.  (38)  can  be  computed  using  the  approximation 

J^(Tr)  s  1  +  (1/4)  (Tr)^  ,  (41) 


when  |T|  r  is  still  smaller  than  unity.  A  correction  factor,  due  to  radial  variations,  is 


pR  ~  p2 

C=j  |j^(Tr)rrdr/^ 

where  for  |  T  |  r  <  1, 


(42) 


C 


+  higher  order  terms 


(43) 


For  all  microwave  frequencies  one  finds,  for  the  case  corresponding  to  Fig.  2,  enormous  correc¬ 
tions,  even  for  a  diode  with  an  unusually  small  radius  such  as  7  |i  which  has  been  used  in  some 
experiments.  Further  numerical  calculations  reinforce  the  preceding  observations  that  the  con¬ 
ventional  approximate  result  for  power  and  impedance,  Eq.  (38),  is  not  adequate  for  the  unysym- 
metric  bulk  modes. 

An  accurate  calculation  of  device  impedance  which  takes  into  account  the  finite  device  radius 
can  be  based  on  the  accurate  calculation  of  power  from  Eq.  (30).  This  will  be  discussed  in  a 
later  work. 
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